Partially suppressed long-range order in the Bose-Einstein condensation of polaritons 
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We adopt a kinetic theory of polariton non-equilibrium Bose-Einstein condensation, to describe 
the formation of off-diagonal long-range order. The theory accounts properly for the dominant role 
of quantum fluctuations in the condensate. In realistic situations with optical excitation at high 
energy, it predicts a significant depletion of the condensate caused by long-wavelength fluctuations. 
As a consequence, the one-body density matrix in space displays a partially suppressed long-range 
order and a pronounced dependence on the finite size of the system. 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) is one of the most 
remarkable manifestations of quantum mechanics at the 
macroscopic scalcji The key feature of BEC of an inter- 
acting Bosc gas is the formation of off-diagonal long range 

the fact that the spatial corre- 
(where n(r, r') is 



order (ODLRO)^ii i.e 
lation g^^\r,T') = n{r, r')/ y/n{r)n{r') 
the one-body density matrix, while n(r) is the particle 
density) extends over the whole system size. The BEC 
mechanism is the key to understand superconductivity 
and superfluidityfi and was the object of renewed inter- 
est following its recent discovery in diluted alkali atomsi^ 
Another candidate system for the observation of BEC 
is that of excitons^ or exciton-polaritons^ in insulating 
crystals. Recently, in particular, many theoretical^iSiiSiiiS 
and experimentaliiii^i'i^'i^ efforts have been devoted to 
BEC of microcavity polaritons. Due to the strong mutual 
polariton interactions, the system is expected to deviate 
significantly from the ideal Bose gas picture. Interac- 
tions are responsible of quantum fluctuations, namely of 
the emergence of collective eigenmodes that differ from 
the single-particle statesi^ Quantum fluctuations in the 
polariton gas are expected to be more relevant than in a 
diluted atomic gas. This will lead to a depletion of the 
condensate that can be very important, reminding of the 
prototypical case of ^Hc in which a condensate fraction of 
less than 10% at equilibrium is achieved. Deviations from 
thermal equilibrium might still enhance this tendency. 
As an example, theoretical predictions for bulk polari- 
tons suggest that an initial condensate population can be 
totally depleted by quantum fluctuations within a suffi- 
ciently long timoi^. In the present full depletion is 
only prevented by the short polariton lifetime, allowing to 
a considerable fraction of the particles in the condensate 
to undergo radiative recombination before scattering to 
the excited states. Another indication of the importance 
of quantum fluctuations comes from the measurement^i 
of the second-order time-correlation function at zero de- 
lay g^^\0). There, quantum fluctuations might explain 
the large deviation from the quantum limit ((/'^'(O) = 1), 
observed far above the condensation threshold. The role 
played by quantum fluctuations is even more important 



in the light of the two-dimensional nature of microcavity 
polaritons. For a two-dimensional system in the ther- 
modynamic limit, as stated by the Hohenbcrg-Mcrmin- 
Wagner theorem^i^ the long-wavelength fluctuations di- 
verge, and ODLRO cannot arise. However condensation 
is prohibited only in infinitely extended systems. Re- 
alistic polariton systems, based both on III-V— and on 
Il-Vli^ semiconductors, exhibit disorder that tends to 
localize the lowest polariton levels over a few tens of 
fim-^. In these situations, it was rigorously proved that 
the discrete energy spectrum, arising either from quan- 
tum confinement in a finite systemic or from disorder 
induced polariton localizatioufi^ allows condensation in 
a thermal equilibrium situation, although the occurrence 
of ODLRO could be inhibited, depending on the local- 
ization length and on the disorder amplitudeF^^ Very re- 
cently, a direct measurement of the spatial correlation 
function in a II- VI semiconductor based microcavity^^, 
provided a striking experimental signature of polariton 
condensation with formation of ODLRO. 

In order to experimentally assess polariton BEC, it is 
important to predict how ODLRO manifests itself in a 
non-equilibrium situation, in a localized geometry, and 
in presence of quantum fluctuations. To this purpose, 
a ficld-thcorctical approach is required for a proper de- 
scription of the relaxation kinetics and of the quantum 
fluctuations in presence of mutual interaction. 

In this work we develop a kinetic theory of polaritons 
subject to mutual interaction, in which the field dynamics 
of collective excitations is treated self-consistcntly along 
with the condensation kinetics. We start from a number- 
conserving Bogolubov approach2£i2i that describes the 
collective modes of a Bose gas properly accounting for 
the number of particles in the condensate. This is re- 
quired in order to develop kinetic equations for the de- 
scription of condensate formation. We then derive a hi- 
erarchy of density matrix equations, including polariton- 
phonon scattering via deformation-potential interaction 
and cxciton-exciton scattering in the exciton-like part 
of the lowcr-polariton branch. This latter mechanism is 
based on the model recently developed by Porras et al^. 
The hierarchy is truncated to include coupled equations 
for the populations in the lower polariton branch and 
for the two-particle correlations between the condensate 
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and the excitations. For the truncation, we assume that 
higher-order correlations evolve much faster than the re- 
laxation dynamics. The kinetic equations obtained in 
this way are solved numerically, assuming a steady-state 
pump at high energy within the exciton-like part of the 
polariton branch. The solution is carried out by account- 
ing sclf-consistcntly for the density-dependent Bogolubov 
spectrum of the collective excitations of the polariton 
gas. We show how this model predicts a dominant effect 
of quantum fluctuations that result in a significant con- 
densate depletion under typical excitation conditions. In 
particular, we discuss the role of quantum confinement 
in a system of finite size and show how ODLRO mani- 
fest itself in typical experimental conditions. Our results 
provide a clear explanation of the partial suppression of 
ODLRO that characterizes the experimental findings^. 

The paper is organized as follows. In Section II we de- 
scribe in detail the theoretical framework of the present 
analysis. Section III is devoted to the presentation of the 
numerical solution of the kinetic equations and to the dis- 
cussion of the results in the light of recent polariton BEG 
experiments. In Section IV we present our conclusions. 



II. THEORY 

We consider the polariton in the lower branch of the 
dispersion as a quasi-particlc in two dimensions)^ de- 
scribed by the Bose field pk, obeying [j5fe,p^,] = 5kk' ■ 
The lower polariton Hamiltonian in presence of Coulomb 
and polariton-phonon scattering is2i 

kk' q 

~\9it{h\ + h^MPk' ^viPk), (1) 

where Hq — TuokP^Pk + ^^q^\^q the free Hamil- 
tonian for polaritons and phonons. The quantity w^^^, 
arises from the Coulomb interaction between excitons 
vxx and from the oscillator strength saturation due to 
Pauli exclusion Wg^j^i^iSl, as 

,(«) 



kk' q 



•^kk' 



VXxXk+qXk'^qXkXk' + 
+VsatXk'-q{Ck+qXk + Xk+qCk)Xk' ■ (2) 

Here the Xk , Ck arc the Hopficld coefficients represent- 
ing the excitonic and the photonic weights respectively, 
in the lower polariton field, i.e. pk ~ XkBk + CkAk- Here 
Bk and Ak are the exciton and the photon destruction 
operators respectively^!. In Eq.Q we consider the in- 
teraction matrix elements in the small-momentum limit 
(consistently with the restriction of the analysis to the 
lower polariton branch of the dispersion)^i22i2^, i.e. 



Vxx 



Vsat 



A 

UsatA' 



(3) 



where Ei, is the exciton binding energy, oq is the exci- 
ton Bohr radius, 2?if2^ is the microcavity vacuum-field 
Rabi splitting, Ugat = 7/167rag is the exciton saturation 
densit y^^i^^ and A is the system quantization area. The 
quantity g^^^, describes the deformation potential inter- 
action of polaritons with acoustic phonons^iS^ and reads 



XkXk' X 



y 2pALzU 

x[ajH\q\)I^{q.)-a^lll{\q\)lt{q.)]. (4) 



Here, p and u are the density and the longitudinal sound 
velocity in the semiconductor, respectively, Lz is the 
quantum well width and ae(fc) are the electron(hole) de- 
formation potentials. The terms and I^(^i^){qz) 
are the superposition integrals of the exciton envelope 
function with the phonon wave function (plane wave) in 
the in-plane and z-directions, respectively, and read 

Ieih)i\l\) = [l + (ao'Z.me(M/2(me-f m;,))2]-3/2 



^ewfe)= / rf^l/ew(^)Pe»'^^ (5) 



where me{h) is the electron (hole) mass while fe{h) (z) 
are the electron(hole) envelope functions in the growth 
direction, according to the exciton envelope function pic- 
turo2£. 

For a kinetic description, a number-conserving 
approach22i2i is required. This formalism allows treating 
in a self-consistent way the field dynamics and the pop- 
ulation kinetics, and describes correctly the condensate 
in- and out-scattering rates. In the number-conserving 
approach, the polariton field is expressed as 



Pk = PkCl+Pk, 



(6) 



i.e. the sum of a condensate operator PkO, and a single- 
particle excitation operator pk<^ The condensate oper- 
ator obeys Bose commutation rules [a, a)] = 1. The 



quantity Nc 



TL^a) is the population of condensed 



particles, while Pk represents the condensate wave func- 
tion in momentum space. The single-particle excitation 
field Pk is orthogonal to the wave function of the con- 
densate, i.e. X]fe PkPk = 0, and obeys the modified 
Bose commutation rule [pk,Pki] = Skk' — PkPk'- Using 
these definitions, the total population at momentum k 
is Nk = {plpk) = \Pk\^N^ + Nk, where Nk - {plpk) is 
the population of non-condensed particles. Within the 
number conserving approach, a quantum fluctuation is 
defined by the operator— 



A 



1 



(7) 



that promotes a particle from the condensate to the ex- 
cited states. N is the total number of particles. This 
field can be formally written via a Bogolubov transfor- 
mation as Afc = UkAk + V*ua^_u, where Uk and V*u are 
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modal functions, and dfe are Bose operator describing the 
collective excitations at energy E/^A^. 

Our model is based on two key-assumptions. First, we 
separate the single-particle energy spectrum into a lower 
energy coherent part and an incoherent part at higher 
energies. This is depicted in Fig. [T] (a) for the typi- 
cal energy-momentum dispersion of the lower polariton 
branch. In a condensate, collective Bogolubov excita- 
tions affect mostly the lower-energy part of the spectrum. 
As the energy becomes larger than the total interaction 
energy vNc (where w is a measure of the interaction ma- 
trix element), Uk 1 and 14- 0, and the collective 
modes thus approach the single-particle states. Hence, it 
is mostly in the low-energy region that quantum fluctua- 
tions will affect the condensate kinetics. A customary^ 
approximation consists in restricting the quantum kinetic 
treatment to the coherent region, while the dynamics 
within the incoherent region is modeled in terms of a 
simple Boltzmann population kinetics. 

The separation, for the lower polariton branch, nat- 
urally coincides with that between the strong-coupling 
region and the flat cxciton-like part of the dispersion, 
as illustrated in Fig. [Ha). The second approxima- 
tion is made possible by the remark that, given the 
large Coulomb scattering amplitude, the field dynam- 
ics of collective Bogolubov excitations takes place much 
faster than energy relaxation mechanisms, made slow by 
the steep polariton energy-momentum dispersion that re- 
duces the space of final states available for scattering pro- 
cesses. We thus assume that, on the timescale of the 
relaxation, a quasi-stationary spectrum of collective Bo- 
golubov excitations arises, that evolves adiabatically and 
is computed self-consistently at each time step in the ki- 
netics, by means of the Popov version of the Hartree- 
Fock-Bogolubov approximation (HFB Popov). This ap- 
proximation is justified by the two following arguments. 
First, in equilibrium conditions, it is known that the HFB 
Popov approximation predicts collective excitations in 
excellent agreement with the measured excitation spectra 
of a weakly interacting Bose gas, when the temperature 
is sufficiently low, i.e. T < 0.5Tc, where Tc is the critical 
temperature. Alternatively, this condition corresponds 
to a density considerably above the critical densityj ^^i'^^ 
Second, our present purpose is not to determine the ex- 
act spectrum of the collective multipole oscillations of 
the polariton gas close to the condensation threshold. It 
is rather to estimate how the density-dependent changes 
in the spectrum affect the relaxation dynamics and the 
coherent scattering processes significantly above the con- 
densation threshold, when ODLRO becomes detectable. 
In addition, close to threshold, the polariton relaxation 
dynamics and the coherent scattering processes are ex- 
pected to be only marginally affected by the the details of 
the spectrum, because the populations in the condensate 
and in the low-lying excited states are small, as argued 
in Ref. m. 

The kinetics is described in terms of a density-matrix 
hierarchy whose time-evolution is obtained from the 



Heisenberg equations of motion. Coulomb interaction 
terms within the incoherent region are treated consis- 
tently with the Boltzmann picture, as was done by Porras 
et al^. In particular, we assume that, in the incoherent 
region, the population of the exciton-like polaritons is 
thermally distributed with an effective temperature de- 
fined by fc^Tj; = ex/ Tlx, where ks is the Boltzmann con- 
stant, Ux is the total particle density and is the total 
energy density in the incoherent region (as shown below, 
these two quantities are determined self-consistently dur- 
ing relaxation, via Boltzmann equations). Within this 
picture, the processes of two particles in the incoherent 
region scattering into one particle in the incoherent re- 
gion and one in the coherent region, result in an effective 
energy relaxation mechanism towards the bottom of the 
polariton branch. In the coherent region, on the other 




wave vector {^m ') 



FIG. 1: (a) Energy-momentum lower-polariton disper- 
sion. Notice the logarithmic horizontal scale, (b) Energy- 
momentum plot of the discrete lower-polariton states in the 
coherent region, as used in the simulations for A — 100 ^rr? . 

hand, we extend the hierarchy of equations to the next 
order, thus including two-body (four operator) correla- 
tions. The most important of these two-body correla- 
tions entering our equations describes coherent processes 
where two polaritons are scattered between the conden- 
sate and the excited states within the coherent region. 
It is given by rhk ~ iV(AfcA_fe) and describes the main 
effect of quantum fluctuations. These processes do not 
conserve energy and could not be described in terms of 
Boltzmann equations for the populations of the single- 
particle states. They are made possible only because, 
in a condensed system, the actual cigcnstates {\v)} are 
collective modes and differ from the single-particle states 
{|fc)}. As a consequence, the quantities {v\KkK-k\u) are 
finite. As already pointed out, their amplitude is ex- 
pected to decrease for increasing energy Ek. We have 
checked that our model reproduces this behaviour, as is 
reported below in connection with Fig. [4] This is proves 
that our approximation, consisting in the separation into 
two energy regions, is consistent with the real behaviour 
of the system. A full treatment in terms of quantum 
kinetics is made here prohibitive by the inclusion of a 
very large number of states in the exciton-likc part of 
the dispersion, that is needed in order to model the ini- 
tial high-energy excitation used in photoluminesccncc ex- 



4 



periments. Models assuming an initial population in the 
lower-energy part of the spectrunti^, can instead treat the 
full spectrum consistently. In fact, in that case, the only 
relevant coherent scattering processes involve states with 
energy smaller than the interaction energy vN (where 
V is the interaction matrix element and N is the initial 
condensate population, see e.g. Fig. 3 of Ref . fisl) . 

The polariton-phonon scattering is treated within a 
shifted-pole Markov approximation, resulting in standard 
Boltzmann contributions^^. For the calculations, we as- 
sume a finite-size homogeneous system of square shape 
and area A with periodic boundary conditions, result- 
ing in spatially a uniform condensate wave function, i.e. 



Jfe.O 
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In a realistic condensate, this assumption 



is valid everywhere, except within a distance from the 
boundary equal to the healing length £, — Ti/ %/ MvriX. For 
polaritons, we find ^ w 1 /im for the estimated density 
threshold. The confinement can model both finite size 
polariton traps^ and the situation close to a local min- 
imum of the disorder potential in extended systemsJ^iii 
The finite size results in quantum confinement and thus 
in a discrete energy spectrum, with a gap between ground 
and first excited state A = {2'it)'^ / {MpoiA)^ 

Within these prescriptions (see Appendix) the kinetic 
equations arc: 



N, = 



rrik 



-7oN, + + N^\^^ + -Y^ 4'l,Im{m,} 



70 + ^^k + ^vl^liNc -Nk~ 5/2) 



m-k 



L q 



2^(1 + 27V,) ^«(^l,(ptpt 



qPkP-k) 



IxTix+n^ \ph + nx \xx + f 



(1 + 2Nk) 



(8) 



The 7fc = 7c|C'a:P is the polariton radiative lifetime, 
|Cfcp being the photon fraction in the polariton state, 
7c the cavity photon lifetime, and the exciton life- 
time. The quantity / denotes the pump intensity produc- 
ing a population in the incoherent region. The suffixes 
"p/i" and "XX" denote the Boltzmann scattering terms 
for polariton-phonon scattering2^ and for exciton-exciton 
scattering in the incoherent region^^. 

In the equation for ifik some simplifications were in- 
troduced. First, wc have fully neglected the contribu- 
tions due to both energy relaxation mechanisms. This 
is consistent with our adiabatic assumption, as the dy- 
namics of the correlation ihk reflects the time-evolution 
of collective excitations, taking place on a much faster 
timescale than relaxation. Second, the higher order 
three-body correlations, arising as the next level of the 
correlation hierarchy, have been factored into products 
of one- and two-body correlations. Furthermore, we 



assume the identity {a)a)aa) = Nc{Nc — 1) to hold, 
as expected for a macroscopic condensate occupation, 
Nc ^ 1. Third, always according to the adiabatic ap- 
proximation, the two-body correlation function between 
condensate excitations, {pljP^_^pkP-k) is evaluated in a 
quasi-stationary limit. It can thus be expressed in terms 
of the modal functions Uk and V*^^ obtained by di- 
agonalizing the stationary Bogolubov problem at each 
time step in the kinetics. Starting from the dynam- 
ical equation Afc {ujk + v'^fiik)^k + vi^l^k^^, with 
£_k = {Nc ~ Nk)/N, we derive the actual eigenvalues 



Ek = [{uJk + v^l^k)' 
functions are given by 



\Vk\=ik- 



"k,k 



while the modal 



Ek-iiVk+vl^l^k) 



Ek-icJk + v):>Ck) 



(0), 



(9) 



and the normalization |C/fep — |V/jp — £,k- In this limit (see 
Appendix) wc can replace in ([5]) the following expression 



{plpl,PkP-k) ^ r{N) 



J2UkVC{l + 2N,) 



+ ^XkNk [xkNk + 2 \Vk\^) + 2 \Vk 

k 



(10) 



where T(iV) 

|2 



N^[{N, + 1){N, + 2)]-i and Xk = Ck + 
2\Vk\^- Nk = {a\ak) is the population of the Bogolubov 
modes, which can be expressed in terms of the singlc- 
particlc population via the exact relation Nk = {N/{Nc + 
1)) X [i\Uk? + \Vk\^)Nk + |14P]. This finally brings to 
a closed set of kinetic equations for the amplitudes m^, 
the populations Nk, and the total density in the 
incoherent region. 

Before discussing the numerical solution of Eqs. (O, 
we present the detailed expressions of the Boltzmann 
terms entering these equations. Their expressions are 
formally identical to the corresponding terms derived in 
Ref. m (for polariton-phonon scattering) and in Ref. [13 
(for exciton-exciton scattering), with an important dif- 
ference: in the present case, the actual HFB Popov spec- 
trum Ek replaces the non-interacting single-particle spec- 
trum. Here we report the expressions of these terms, the 
details of the derivation being described in RefsiSii^. 



A 



W: 



ph 



^nxil + Nk)-W^txNk{v+'-f) 



Nk 



E [w!tkNk'{l + Nk) + 

k'eUaah 

-WPtk'Nka + Nk')] , 



Nk 



XX 



W^J^k^lil + Nk) - W^J^^n^Nk, 



(11) 



(12) 
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(13) 



and 



\xx 



A ^ 

keUaoh 



(14) 

where ijA = [(mg + 'mii)kBTx]A/2TTfi is the number of 
states in the incoherent region, and Ucoh denotes the co- 
herent region. Here, Nk = Nk for fc ^ and A^o = ^c- 
Note that in Eqs. (|imi4p we use the fact that n^j/ry <C 1 
and k <^ kx, where kx is the averaged momentum of the 
exciton distribution in the incoherent region. The rates 
appearing in Eq. (jlHll4p are given by 



W, 



ph 



k^k' 



h 



9k,k' 



^ \Ek - Ek 



ph 
^kk' 



B(Ek-Ek') 



= ^\k^W+^z = l^^fe - Ek' \lht 



(15) 



where is the population of phonons having energy 



E = Eh 



Eh 



AH 



' ji2^(Ek-Ek^)/kBT^ 



w. 



XX 



[rrie 



■mil. 



"^A, A, 



,2{E^-E^^)/knT^ 



(16) 



(17) 



In all these expressions, consistently with our assump- 
tion the spectrum in the incoherent region Ek^ is the 
single-particle spectrum, only accounting for the overall 
density-induced blue shift of the polariton branch. Be- 
sides the set of equations © , the model is completed by 
the introduction of an equation describing the evolution 
of the energy density in the incoherent region Cx- Fol- 
lowing the procedure of Ref. [2^ . this equation reads 



EEk r 
~A 



keu^_oh 

1. (kBTx) Ux + {kBW f - VJP'' 



(18) 



Here, the first term represents the heating of the inco- 
herent population, produced by the exciton-exciton scat- 
tering process and imposed by energy conservation; the 
second term is the cooling due to the exciton radiative 
recombination; the third term is originated by the as- 
sumption that the incoherent population is created at 
the lattice temperature T^; the fourth term represents 
the cooling induced by exciton-phonon coupling and it is 
evaluated as in Eq. (21) of Ref. [H. 



III. NUMERICAL RESULTS 

Numerical solutions assuming a steady state pump 
have been computed in the time domain. We assume 



parameter values relative to the CdTe microcavity of 
Ref. [13, with Rabi splitting 2hQB. = 26 meV and cav- 
ity photon-exciton detuning 5 = 5 meV, at the lattice 
temperature T = 10 K. The quantization area is as- 
sumed everywhere A = 100 /.jm^, unless specified, con- 
sistent with estimates of polariton localization lengthfi^ 
and gives rise to the discrete set of polariton states plot- 
ted in Fig. [ijb). For this system we obtain from Eq. 
dg) vxx = 3.3 X 10-5 meV, v,at = -0.5 x 10-'^ meV 
and the resulting polariton-polariton interaction matrix 

element at zero momentum is Wg",] = 6 x 10^^ meV. The 
solutions always display steady-state long-time values af- 
ter an initial transient, as shown in Fig. Here we 
also notice that, during the early stages of the conden- 
sate growth, the scattering processes favor condensation, 
through the positive values taken by Imjmfe}. immedi- 
ately afterwards, when a large condensate population is 
reached, the quantity Im{mfe} changes sign and coher- 
ent scattering terms start depleting the condensate. In 
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400 



FIG. 2: Time dependent results for the populations of the 
condensate state Nc and of the first excited state A^'i and for 
the imaginary part of the coherent scattering amplitude mi, 
for a given value of the excitation pump / above condensation 
threshold. After an initial transient, all the quantities reach 
a stationary value. 

Fig. [3] (a), we plot the steady-state populations per state, 
for varying pump intensity. A pump threshold is found 
at about / = 12 ps"^ /im"^, corresponding to a polari- 
ton density n = N/A ~ 10 ^m~^ and a total exciton 
density Ux — 100 /im^^ in the incoherent region. Let 
us remind that the system studied in Ref. [ij is com- 
posed by 16 quantum wells. As the polariton modes are 
expected to have similar amplitudes at all the quantum 
well positions, in the experimental literature the exciton 
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density per quantum well is usually estimated by sim- 
ply dividing the total exciton density by the number of 
quantum well o^^i^^ . In the present case, following this 
prescription, the exciton density per quantum well is 7 
/im~^, i.e. two orders of magnitude lower than the ex- 
citon saturation density in CdTe risat — 5 x 10'^ iimr'^ . 
This confirms the validity of our description of polaritons 
in terms of a weakly interacting Bose field. 

Above threshold the condensate population becomes 
macroscopic. Its growth for increasing / is however sup- 
pressed by the corresponding increase of the population 
of low energy excitations. Consequently, the population 
distribution cannot be fitted by a Bosc-Einstcin func- 
tion. The discrepancy is partly due to the Bogolubov 
quasiparticle spectrum - characterizing an interacting 
Bose gas at thermal equilibrium - and partly to quan- 
tum fluctuations. In order to distinguish the two con- 
tributions, we compare the kinetic result to a distribu- 
tion computed for an equilibrium interacting Bose gas 
in the HFB Popov limit accounting for spatial con- 
finement. For this equilibrium distribution, we assume 
the same density as obtained from the kinetic model for 
a given pump /, while the temperature is extrapolated 
from the slope of the high-energy tail of the same ki- 
netic model distribution. In Fig. [3] (a), the result for 
/ = 50 ps^^ /_im^^, is compared to the equilibrium HFB 
Popov distribution with n = 100 /im~^ and T = 20K. 
As expected, the equilibrium result already deviates from 
the ideal distribution, due to the modified spectrum of 
the interacting system. However, equilibrium and kinetic 
results differ significantly in the low-energy region. In 
particular, the kinetic model predicts a larger conden- 
sate depletion. The difference is due to the dominant role 
played by quantum fluctuations (see also Fig. ([4]) below), 
whose amplitude deviates from the equilibrium predic- 
tion and has to be evaluated by means of a proper kinetic 
treatment like the present one. Also in the kinetic model, 
the energy dispersion is modified by the presence of the 
condensate because of the two-body interaction, display- 
ing the linear Bogolubov spectrum of collective excita- 
tions at low momenta,— as shown in Fig. [3] (b). Here, 
the dashed line highlights the linear part of the spec- 
trum. The plot shows that, even at the highest pump 
intensity, this feature extends over an energy interval 
smaller than 0.5 meV, well within the measured spectral 
linewidthpi^ Samples with a significantly smaller polari- 
ton linewidth (longer radiative lifetime) would therefore 
be needed, in order to measure this distinctive feature 
of BEC. In Fig. U] we show the values of the coherent 

scattering rates v^^\fjjm{'mk\ as a function of the energy 
of the corresponding states. As expected, they decrease 
dramatically for increasing energy and their contribution 
to the dynamics vanishes for states outside the coherent 
region, thus confirming our initial assumption of separa- 
tion into two momentum regions. We see also that, for 
increasing excitation intensity the values of the coherent 
scattering terms increase, resulting in an increased am- 
plitude of quantum fluctuations. 
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FIG. 3: (a) Steady state populations for increasing pump 
intensity /. Open squares: equilibrium HFB Popov solution 
for (n,T) corresponding to the steady-state solution for / = 
/a. Thin line: B-E distribution fitted to the high-energy tail 
and to the condensate population for f = fa- (b) Energy 
dispersion below and above threshold (same legend as above) . 
The dashed line is a guide to the eye to highlight the linear 
part of the dispersion. 
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FIG. 4: Coherent scattering contributions ii^'^^^Im{mfe}, for 
two values of the excitation pump /. The amplitude of such 
processes decreases dramatically for increasing energy and so 
their contribution to the dynamics vanishes outside the co- 
herent region. 



For increasing system area A, the condensate frac- 
tion in the steady-state regime decreases, as shown in 
Fig. [5] (a), because coherent scattering is favored by a 
smaller energy gap A. Thermal and quantum fluctu- 
ations will eventually dominate in the thermodynamic 
limit of infinite size, resulting in a full condensate de- 
pletion, as required by the Hohenbcrg-Mcrmin- Wagner 
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FIG. 5: (a) Condensate fraction as a function of the system 
area for / = 30 ps~^ /xm~^. (b) First order spatial correlation 
function below and above threshold, compared to the same 
quantity resulting from the equilibrium HFB Popov model. 



in an experiment in which the light emitted by differ- 
ent positions on the sample is made interfere. For such 
an experiment, we therefore predict the increase of the 
spatial correlation length as a signature of condensation, 
with a correlation staying below 0.5 because of quantum 
fluctuations, even far above the threshold. 

Finally, we mention that we have also performed sim- 
ulations for material parameters modeling a GaAs mi- 
crocavity with Rabi splitting 2hD,ji ~ 7 meV, detuning 
(5 = meV, lattice temperature T = 10 K, and system 
area A = 100 iim^. In this case vxx = 6 x 10~^ meV, 
Vsat = —0.15 X 10~^ meV and the resulting polariton- 
polariton interaction matrix element at zero momentum 



IS v. 



(0) 
0.0 



1.5 X 10^^ mcV. For this case as well, we ob- 
serve the occurrence of polariton condensation and the 
partial suppression of the ODLRO (not shown) . Quanti- 
tatively, we notice that in this case the total exciton den- 
sity at threshold is rix — 500 /im~^, significantly higher 
than in the previous case. Nevertheless, as the number 
of quantum wells needed to achieve 7 meV Rabi splitting 
is at least 4, the estimated exciton density per quan- 
tum well is about 100 /im~^, namely still one order of 
magnitude lower than the saturation density for GaAs, 
Usat = 2 X 10'^ /xm~^i^ Our approach is thus justified 
also in this case. 



IV. CONCLUSIONS 



theorem. Polariton condensation occurs thanks to the 
locally discrete nature of the energy spectrum, induced 
either by artificial confinement or by disorder. In a real- 
istic systeroi^iiii^^, localization could therefore affect the 
polariton BEG, independently of other parameters like 
Rabi splitting and exciton saturation density. 

We finally study the influence of quantum fluctua- 
tions on the one-body spatial correlation ^'^^(r, r') = 
n{r, r') / •\/n(r)n(r'). The one-body density matrix 
n(r,r') is the direct expression of ODLRO that char- 
acterizes BEG.i^ It depends on |r — r'| for a uniform 
system and can be computed in terms of the Fourier 
transform of the population Nk- The density in the 
denominator renormalizc the shape of the condensate 
wave function, hence we expect the averaged quantity 
5(1) (r) = l/A jd'Rg^^^R,'R + r) to be scarcely affected 
by the assumption of a uniform condensate. In Fig.[5](b), 
we plot g^-^^ (r) below and above the condensation thresh- 
old. Below threshold, correlations vanish for distances 
larger than 1 — 2 fim, as predicted by both the kinetic 
model and the equilibrium HFB Popov approach. Above 
threshold, the correlation length increases and g^^\r) re- 
mains finite over the whole system size. However, for all 
values of the pump, ^^^■'(r) remains smaller than 0.5 at 
distances larger than 5 /xm, whereas the equilibrium HFB 
Popov result is significantly larger, due to a larger con- 
densate fraction. Therefore the effect of quantum fluc- 
tuations turns out to strongly affect the formation of 
ODLRO. The quantity g^^\r) could be easily accessed 



In conclusion, the present model shows that the dy- 
namics of quantum fluctuations significantly affects po- 
lariton BEG and the formation of ODLRO in a polari- 
ton condensate. In a typical case, quantum fluctuations 
partially deplete the condensate, already slightly above 
threshold. Quantitatively, the effect depends on the lo- 
cally discrete energy spectrum, due to trapping or to 
disorder. We predict that the observation of BEG and 
ODLRO should be favored by smaller polariton size, as 
in the recently studied polariton "quantum boxes" ^ or 
in local minima of the disorder potential. This suggests 
that, for a given sample, a study of the polariton local- 
ization length in the lowest energy statesii could give 
deeper insight into the BEG mechanism. 
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APPENDIX 

Equations ([5]) arc derived under the assumption that 
3-body and higher order correlation terms can be fac- 
tored. In this Appendix we provide some detail on their 
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derivation. Equations for the populations Nc and Nk, 
with respect to the polariton-polariton interaction part, 
result directly from Heisenberg time-evolution of oper- 
ators. The corresponding exciton-phonon term N\ph is 
obtained within a standard Boltzmann kinetics in the 
Markov limit ^, while for the exciton-exciton term N\xx 
we have used the equations introduced by Porras et al. 
in Ref. H^- Some less straightforward steps were taken 
to obtain the equations describing the scattering ampli- 
tudes TTife. First, as pointed out in the text, we included 
only the dynamics induced by polariton-polariton scat- 
tering Hamiltonian within the coherent region. Then, 
Heisenberg equations of motion result in the following 
equations 

ihrhk = 2{hujk - v'i'l)rhk + v'^l{a^ a\a^ a + aa))pkP-k) 



From here, within the assumption of a spatially homoge- 
neous system, som of the momentum sums can be carried 
out explicitely and Eq.([8]) is finally obtained. 

We now turn to the evaluation of the two-body corre- 
lations (pJpLqPfep_fc). As explained in the text, they are 
evaluated in a quasi-stationary limit, via the stationary 
solution of the Bogolubov problem. In terms of collective 
excitation operators , defined in Eq. ([7]) , this quantity 
is rewritten as: 



iPlp'LgPkP-k) 



[aaa)a)p\p^_qPkP-k) 



{N, + l){N, 



(iVe+l)(iVc + 2) 



2) 

(AtAL,AfcA_,.). 



(A.7) 



+ J2''g'^^''''^''''"-''^Pl+q'-kP-k+PkPl+q'+k)PlPq') 
qq' 



The correlation amplitude for the field A^ is then com- 
puted by using the Bogolubov transformation Afc = 



Ukak + V*a[ 



All the resulting terms are factored 



as products of the Bogolubov quasi-particlc populations 



Nk 



c^lctk), 



Then, by neglecting 3-body correlations, we can intro- 
duce the following factorizations: 



{al&lakaq) ~ Nk {Nq - Skq) ■ (A.. 



{a) a) a) apkP-k) ^ (A^c - 2)mfe; 



(A.2) 



{a}a^p\pqPkP-k) ^ {Nq - 5q^kNk - Sq,^kN-k)mk 



+Nq^kmq-k + Nq-kmq.k] 

{a^ap\p^_qPkP-k) ^ Nc(plplqPkP-k); 
{a'<a^aaplpk) N,{N, - l)Nk; 

{a^a^pl+q,-kP-'kPqPq') ^ 2N q+q, ^k^q' rhq^^k 
+ Nq^qf^k,-krh'q,q> 
-Sq.kSq' .kNkrhk). 



(A.3) 
(A.4) 
(A.5) 



(A.6) 



Collecting all the terms, wc finally recover Eg. pT]]) . In 
this equation, single-particle and quasi-particle popula- 
tions are related with each other through the expression 



(plpk) ^ j^iiw + m'Wk + \Vkn (A.9) 



N 
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